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ABSTRACT 



We have obtained optical long-slit spectroscopy of the nucleus of M32 
using the Space Telescope Imaging Spectrograph aboard the Hubble Space 
Telescope. The stellar rotation velocity and velocity dispersion, as well as the 
full line-of-sight velocity distribution (LOSVD), were determined as a function of 
position along the slit using two independent spectral deconvolution algorithms. 
We see three clear kinematical signatures of the nuclear black hole: a sudden 
upturn, at ~ 0.3" from the center, in the stellar velocity dispersions; a flat or 
rising rotation curve into the center; and strong, non-Gaussian wings on the 
central LOSVD. The central velocity dispersion is ~ 130 km s" 1 (Gaussian 
fit) or > 175 km s^ 1 (corrected for the wings). Both the velocity dispersion 
spike and the shape of the central LOSVD are consistent with the presence of a 
supermassive compact object in M32 with a mass in the range 2 — 5 x 1O 6 A^0. 
These data are a significant improvement on previous stellar kinematical data, 
making M32 the first galaxy for which the imprint of the black hole's gravitation 
on the stellar velocities has been observed with a resolution comparable to that 
of gas-dynamical studies. 

Keywords: galaxies: elliptical and lenticular — galaxies: structure — galaxies: nuclei 
— stellar dynamics 
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1. Introduction 

The presence of a supermassive compact object, presumably a black hole, at the center 
of the dwarf elliptical galaxy M32 has been suspected for some time ( |l'onry 1987] ). The 
evidence consists of rapid rotation of the stars near the center of M32 and a central peak in 
the stellar velocity dispersions flTonry 19871 ; pressler fc Richstone 19881 ; ICarter fc Jenkins] 



|1993| ; [van der Marel et al. 1994a] ; Pender, Kormendy fc Dehnen 19961) ■ The most recent 



study (|van der Marel et al. 1997] , 1998) used data from the Faint Object Spectrograph 
(FOS) on the Hubble Space Telescope (HST) to infer the rotation and dispersion velocities 
with a spatial resolution of ~ 0.1" in the inner ~ 0.5" of M32. The FOS data revealed a 
sharper rise in the stellar velocity dispersions than had been observed from the ground; 
however the velocity resolution of the FOS is limited, making that instrument only 
marginally useful for the study of a low velocity dispersion system like M32. In fact, van 
der Marel et al. (1998) found large point-to-point variations in their velocity dispersion 
measurements, making the dynamical interpretation uncertain. Nevertheless, the case for a 
supermassive black hole in M32, of mass Mh ~ 3 x 1O 6 A^0, was considerably strengthened. 

Here we present observations of M32 carried out using STIS, the Space Telescope 
Imaging Spectrograph, on HST. Our data were obtained as part of the STIS Investigation 
Definition Team's (IDT) key program to observe the nuclei of a sample of ~ 15 nearby 
galaxies in the spectral region centered on the Calcium triplet, A ~ 8600A. This paper is 
the first in a series that will present stellar-kinematical evidence from STIS for the presence 
(or absence) of supermassive black holes in galactic nuclei. 

The STIS data improve on earlier ground-based and FOS data from M32 in several 
ways. The spatial resolution of STIS is ~ 0.1", or ~ 0.3 pc at the distance of M32, similar 
to that of the single FOS aperture; however STIS provides continuous spatial sampling 
along a slit. The spectral resolution of STIS in the G750M mode (~ 38 km/s) is much 
greater than that of the FOS making STIS a more suitable instrument for observing M32, 
whose velocity dispersion outside of the nucleus is only ~ 60 km s _1 . We were able to 
obtain from the STIS spectra not only the lowest-order moments of the stellar velocity 
distribution - the rotation velocity and velocity dispersion - but also the full line-of-sight 
velocity distribution (LOSVD) as a function of position along the major axis. 



A number of modelling studies ( [van der Marel et al. 1994b| pehnen 1995| ; |Qian et 



al. 1995Q have made predictions about the observable signatures at HST resolution of a 
supermassive black hole in M32. The black hole is expected to be associated with three 
kinematical features. (1) The stellar rotation velocity should remain flat, or rise slightly, 
into the central resolution element. (2) The stellar velocity dispersion should exhibit a 
sudden upturn at a distance of ~ 0.2" — 0.5" from the center, reaching a central value of 
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~ 120 km s _1 or greater depending on the mass of the black hole. (3) The distribution of 
line-of-sight stellar velocities in the resolution element centered on the black hole should be 
strongly non-Gaussian, with extended, high-velocity wings. Part of the predicted rise in 
the velocity dispersion near the center would be due to these wings; part to blending of the 
rotation curve from the two sides of the galaxy; and part to an intrinsic rise in the random 
velocities. 

We see all three signatures of the black hole in the STIS data. The velocity dispersion 
spike is most impressive; in terms of the usual parametrization o"o (the dispersion of the 
Gaussian core of the LOSVD), the central measured value is ~ 130 km s _1 , while the value 
of a corrected for the non-Gaussian wings of the LOSVD is considerably greater, at least 
175 km s _1 , or ~ 3 times the value in the main body of M32. 

A detailed description of the observations and the data reduction is given in §2. In §3 
we describe the methods used to extract the stellar LOSVDs from the STIS spectra. We 
carried out independent analyses based on two spectral deconvolution routines: the FCQ 
algorithm of R. Bender (1990), which is a Fourier method; and the MPL routine of D. 
Merritt (1997), which is based on nonparametric function-estimation techniques. The two 
algorithms gave consistent results for the low-order velocity moments of M32 (the mean 
and dispersion) but there were systematic differences in the recovered values of hi, the 
parameter that measures symmetric deviations of the LOSVD from a Gaussian. We argue 
that the h± values are consistent after taking into account the biases in the two methods; if 
our analysis is correct, the true values of in M32 are significantly greater than zero both 
within the region where the black hole's force is dominant and at larger radii. 

A preliminary interpretation of the M32 kinematical data is presented in §4, where we 
address two basic questions: Have we resolved the central spike in the velocity dispersions? 
and: Is the shape of the central LOSVD consistent with what is expected for a nucleus 
containing a black hole? We answer "yes" to both questions, although the resolution of 
the central spike may be marginal if the black hole mass is near the lower limit of the 
allowed range. In any case, M32 is the first galaxy for which the imprint of the black hole's 
gravitation on the stellar velocities has been observed with a resolution comparable to that 
of the gas-dynamical studies (e.g. |Bower et al. 1998| [Ferrarese fc Ford 1999|) . 

We show that both the velocity dispersion spike and the wings of the central LOSVD 
are consistent with the predictions of simple dynamical models containing black holes, with 
masses M h in the range 2 x 10 6 A^ Q < M h < 5 x 1O 6 A^ . More accurate estimates of 
Mh will be presented in Paper II where the full kinematical data set will be compared to 
three-integral axisymmetric models. 
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We adopt a distance to M32 of 0.7 Mpc; thus 1" corresponds to 3.25 pc. 



- 6- 



2. Observations and Data Reduction 

M32 was observed on 1998 September 04 with STIS in the long-slit mode with 
wavelength centered on the Ca II triplet feature near 856 lA. STIS is described by Woodgate 
et al. (1998) and its on-orbit performance by Kimble et al. (1998). The present data are 
part of a survey of the nuclei of nearby galaxies being conducted by the STIS IDT (HST 
Program ID: 7566). The goal of the survey is to place stellar kinematical constraints on 
the masses of nuclear black holes. Two orbits of data with a total integration time of 
4,898 seconds were obtained. The HST spacecraft tracking was operated in fine lock with 
a reported jitter of no more than 3 mas rms or 11 mas peak-to-peak. The aperture was 
52" x 0.1" with a position angle of 163°, coincident with the M32 isophotal major axis. The 
CCD data were read out in the unbinned mode. Spatial sampling at the focal plane was 
at every 0.05071", corresponding to a 2-pixel optical resolution of about 0.115" FWHM. 
Outside of ±0.7", the data were binned spatially to enhance the signal. The spectral 
resolution was approximately 38 km s^ 1 . M32 spectra were obtained in the CR-split mode 
to assist with cosmic ray (CR) identification and rejection. The location of the galaxy center 
in the image was moved by approximately 4.5 rows along the aperture between the two 
orbits to ensure that residual detector sensitivity variations, that may not be completely 
removed from the data during reduction, are not mistaken for weak features. This form of 
dithering also assists with the identification of hot pixels in the CCD which do not rectify 
well. 

A spectrum of the star HR 7615 was obtained with the same STIS aperture earlier in 
the program (Figure [l]). This bright, K0 III giant was used as a template when inferring 
the LOSVD of the stars in M32 from the STIS spectra. Template spectra of a handful of 
other bright, cool stars of various spectral types were also obtained and used to test the 
sensitivity of the spectral deconvolution to template mismatch. A set of spectral images of 
HR 7615 were taken centered at —0.05", 0.00" and +0.05" with respect to the centerline of 
the aperture (i.e. offset along the dispersion direction). These data were added together 
using appropriate weights to match the aperture illumination profile of M32. Each of 
these spectral images has a slightly different shift in its velocity centroid and the combined 
image provides a more accurate template for determining the kinematics of an extended 
object such as a galaxy. Internal wavelength calibration images ("Wavecals") as well as 
an internal continuum lamp image ("Flat field") were taken for calibration purposes. An 
image that has had all of the instrumental response removed is said be a rectified image. 
The contemporaneous flat field spectrum was obtained in the portion of the orbit where 
HR 7615 was behind the earth. The flat spectrum was taken through the 0.2" x 0.09" 
aperture rather than through the long aperture since the former is superior for removing 
the pixel-to-pixel detector response of a stellar point source. A contemporaneous flat field 
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was used to remove the internal fringing which is significant for wavelengths greater than 
7500 A and which changes over time ( |Goudfrooij , Baum fc Walsh 1997]) . We note that the 
fringing is far more serious for low dispersion spectra especially for high S/N than it is for 
our data taken in medium dispersion. 

Cosmic rays (CRs) account for approximately 20% of the total signal and contaminate 
approximately 5% of the pixels in a typical exposure. CRs were identified and removed using 
the following procedure. The centroids in the cross-dispersion direction were determined for 
each rectified image and the images shifted so that the galaxy core appeared on the same 
row. CRs were identified and removed by comparing the flux in a given pixel to the flux in 
the corresponding pixel in subsequent images. For each pixel, outlying values were rejected 
and excluded when the frames were averaged together. Our data set for M32 included 4 raw 
images. Most pixels were found to have 4 frames contributing to their average values, while 
fewer pixels had 3 or less frames. Only 2 pixels within the central 2" had contributions from 
no frames. Those pixels were assigned values representing the average of their adjacent 
pixels. 

The M32 data were reduced using two separate approaches: 1) "Shift and Add," 
and 2) "Frame by Frame." The latter relies heavily on a standard software package 
called CALSTIS at Goddard Space Flight Center. The "Shift and Add" method starts 
by removing the detector response using contemporaneous flat, bias, and dark calibration 
files. CR hits are then removed using the procedure described previously. The frames are 
averaged together and the resulting frame is remapped to place the spectra from a single 
location along the aperture onto a single row. As with most spectrographs, STIS produces 
spectra with S-shaped and pincushion distortions as well as spectra that are not aligned 
exactly with a row. A cubic interpolation was used to remap the spectra for later analysis. 
The remapping is not perfect, with centroid errors of approximately 0.1 pixel rms. This 
level of accuracy was deemed adequate, although we are working to improve it. In addition, 
the remapping produces very minor (< 1%) residual Moire ripples in the data, which can 
be minimized but not completely eliminated. Further work is also underway to measure 
and correct for this residual affect. The strength of the "Shift and Add" technique is that it 
preserves photometric accuracy. However it has the disadvantage of introducing a subpixel 
image smearing since each frame is registered to the nearest integer pixel. 

The philosophy of the "Frame by Frame" method is to apply all calibration corrections 
(including the remapping described above) to a single frame before the resulting frames 
are added together. This approach has the advantage of preserving the highest spatial 
resolution. However, a substantial fraction of the pixels in each frame in this case are 
interpolated values, potentially sacrificing some photometric precision. As noted above, 
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CRs contaminate approximately 5% of the pixels in a 20 minute exposure. After remapping, 
as many as 20% of the pixels must be assigned a reduced weight for the final averaging 
since a single pixel often gets remapped partially into several pixels in the new image. 

Fortunately, both approaches to the data reduction were found to give very similar 
results. While the "Frame by Frame" method produced a higher apparent S/N ratio than 
the "Shift and Add" method, the former method is more prone to the introduction of 
systematic error. We adopted the more conservative approach of accepting a somewhat 
higher variance rather than risk the introduction of a bias. We therefore adopted the " Shift 
and Add" spectra for the analysis of the LOSVDs in the present study. Figure 1 shows 
spectra at several positions along the aperture. As noted above, we are continuing in our 
efforts to refine the data reduction techniques still further. 
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3. Recovery of the Stellar Velocity Distribution 
3.1. Method 

An observed spectrum /(A) is the convolution of the line-of-sight velocity distribution 
N(V) of the stars within the aperture with the spectrum of a single star T(X): 



/(In A) = J N(V)T(\n\- V/c)dV. 



The goal is to extract an estimate N(V) of the true stellar broadening function N(V) given 
/(A) and T(A), both observed with the same instrument. For T(A) we adopt the template 
spectrum of Figure ffl. 

Two independent deconvolution routines were used. The first algorithm, the "Fourier 
Correlation Quotient" (FCQ) method ( [Bender 19901 ; |Bender, Saglia fc Gerhard 1994|) , 



constructs an estimate of the broadening function using Fourier techniques. The FCQ 
routine differs from earlier Fourier algorithms (e.g. [Sargent et al. 1977| ) in that the 



deconvolution is based on the template-galaxy correlation function rather than on the 
spectra themselves. This approach is less sensitive to template mismatch QBender 1990| ). 



The second algorithm, "Maximum Penalized Likelihood" (MPL), finds N(V) as the solution 
to a penalized likelihood problem. The MPL estimate of N(V) is computed on a grid in 
V in such a way as to optimize the fit of the convolved template to the galaxy spectrum, 
subject to a "penalty" that measures the lack of smoothness of N(V) ( Merritt 1997|) . 



Both algorithms are nonparametric in the sense that no explicit constraints are placed 
on the functional form of N(V). However they differ in two ways that are important for 
the current study. The FCQ algorithm requires that the absorption lines in the template 
spectrum be narrow compared to the broadened lines of the galaxy spectrum, i.e. that 
the galaxy velocity dispersion be large compared to the instrumental resolution. The 
MPL routine works well even when the galaxy velocity dispersion is small, as long as both 
template and galaxy spectra are observed at the same spectral resolution, at least in the 
case that the template star and galaxy have the same intrinsic absorption line properties 
(an assumption that will not be tested here). 

The two algorithms differ also in the way they deal with the amplification of noise 
that accompanies the deconvolution. The FCQ routine uses a Wiener filter to suppress 
high-frequency components of the template-galaxy correlation function Kt,g- The degree of 
smoothing is determined by a factor, called here W, which fixes the width of the Gaussian 
function used to model the low-frequency, or signal, component of Kt,g- The choice W = 1 
corresponds to "optimal" filtering and larger values produce less smoothing; the FCQ 
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algorithm adopts a default value of W — 1 but automatically increases W (to a maximum 
of 1.3) if the recovered LOSVD shows evidence of significantly non-Gaussian wings. In 
the MPL algorithm, the level of smoothing is determined by a factor a that multiplies the 
smoothness penalty function. This penalty function defines any N(V) that is Gaussian as 
"smooth," regardless of its mean or dispersion, via Silverman's (1982) prescription; in the 
limit a —>■ oo, the MPL estimate of N(V) is the Gaussian function which is most consistent, 
in a maximum-likelihood sense, with the galaxy spectrum. There is no a priori way of 
computing the optimum value of a in the MPL algorithm, a point that we return to below. 

The different effects of smoothing on the form of N(V) are illustrated in Figure |2], which 
shows estimates of the LOSVD in the central resolution element of M32 as computed by the 
two routines. Both algorithms produce rapidly-fluctuating solutions when undersmoothed, 
a consequence of the amplification of noise that always accompanies deconvolution. The 
only significant difference in this regime is the non-negativity of the MPL estimates, a 
consequence of the logarithmic form of the penalty function ( [Silverman 1982j ). As the 



smoothing is increased, systematic differences begin to appear which are related to the 
different smoothing algorithms in the two codes. Solutions obtained via MPL tend to be 
more robust with respect to the degree of smoothing, producing in the limit of large a 
a Gaussian fit. However Figure suggests that estimates of certain quantities, e.g. the 
wings of the LOSVD, might depend sensitively on the choice of smoothing level in either 
algorithm. 

Once an estimate of N(V) has been obtained, various quantites related to the 
line-of-sight velocity distribution can be derived. The simplest of these are the mean 
and rms velocities, which we denote by V and a respectively. As is well known, both 
quantities are difficult to estimate for numerically-recovered LOSVD's since they are 
sensitively dependent on the form of N(V) at large velocities where this function is most 
poorly determined. A standard alternative is to describe N(V) by a Gram-Charlier or 
Gauss-Hermite (GH) series, the product of a normalizing Gaussian with a sum of Hermite 



polynomials Hi, both expressed in terms of (V — Vo)/o"o (Thompson fc Tapia 1990|) . The 
parameters Vo and <To take the place of V and a; while their definitions are to an extent 
arbitrary, these parameters are typically determined by requiring the coefficients of Hi and 



if 2; called hi and h 2 , to be zero ( Gerhard 1993 ; van der Marel fc Franx 1993|) . Because 



Vq and o"o describe the Gaussian core of the LOSVD, they are relatively insensitive to 
deviations of N(V) from Gaussianity at high velocities. Information about these deviations 
is contained in the higher-order coefficients h 3 , /i 4 etc.; h 3 measures asymmetries in N(V) 
and measures the strength of symmetric, non-Gaussian wings. 

The FCQ and MPL algorithms derive the GH parameters from N(V) in slightly 
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different ways; details are given in Appendix A. 

When applied to the STIS M32 spectra, the two algorithms were found to give 
consistent results for the lowest moments of N(V), i.e. Vq, ao and h 3 . However the estimates 
of h,4 differed significantly at positions outside of the central ~ 0.1". The FCQ algorithm 
gave —0.15 < < at almost all positions; negative values of imply an N(V) that 
falls off more sharply than a Gaussian at large velocities. The MPL algorithm gave values 
for /i 4 in the range < < 0.1, almost all positive, corresponding to LOSVD's with 
super-Gaussian wings. Positive values of are expected near a black hole ( |Bahcall fc Woll 
|1976| ; |van der Marel 1994j ) and are also characteristic of models with radially-anisotropic 
velocity distributions. 

We discuss the origin of this discrepancy in Appendix B. We believe that the primary 
reason for the systematic difference in values is the low velocity dispersion of M32. When 
a galaxy's velocity dispersion is comparable to the dispersion of the template star spectrum 
(~ 50 km s _1 in the case of HR7615), the FCQ algorithm has difficulty recovering the true 
LOSVD (Figures B2,3; |Bender, Paquet fc INieto 1991| ). The N(V) 7 s recovered by FCQ in 
this regime are more sharply truncated than the true iV(V)'s, leading to systematically low 
estimates of h^. For values of do and S/N comparable to those of M32 at ~ 1", Figure B2 
shows that the estimates of generated by FCQ depend only weakly on the true /14, with 
a bias that approaches —0.1 for a true /14 of 0.1. The MPL algorithm suffers much less from 
this bias (Figures B3-5). 

Bias in nonparametric function estimates can always be reduced by reducing the 
degree of smoothing (e.g. [Silverman 1986| ), which in the case of the FCQ algorithm means 
increasing W. Figure B5c suggests that increasing W from its default value of 1 to values of 
~ 2 can reduce the bias in FCQ estimates of /14 by factors of 2 or greater, even when a is 
as large as 100 km s _1 . We carried out this experiment with the STIS data; the results are 
shown in Figure [|. The average h.4 values recovered by FCQ in M32 are indeed dependent 
on W; a change in W from 1 to 1.5 has the effect of increasing (/i 4 ) from ~ —0.08 to 
~ +0.08. The latter value is essentially identical to the mean value of h 4 recovered via 
MPL. 

We note that the dependence of /14 on W could be due either to the suggested 
explanation, i.e. the need to include more frequency channels for dispersions close to 
the instrumental resolution, or alternatively to a mismatch between galaxy and template 
spectral properties in the wings of the lines which results in incorrect continuum subtraction. 
We will not explore the second possibility here but note again that the Monte Carlo 
experiments in Appendix B suggest that values of W of order 2 are appropriate even when 
there is no template mismatch. 
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The values of /14 recovered by MPL are also dependent on the value of the smoothing 
parameter a but much less so (cf. Figure 0), until a is made so large that the LOSVD is 
forced into a Gaussian shape. The Monte Carlo tests summarized in Figure B5b suggest 
that the bias in as recovered by MPL is likely to be only of order ~ —0.02, several times 
smaller than with FCQ. 

We conclude that the values of /14 recovered by the two algorithms are consistent once 
their different biases are taken into account and that the values returned by MPL are likely 
to be more accurate. Henceforth we adopt the MPL estimates. 

The full set of GH parameters derived from the STIS spectra and their la confidence 
intervals are given in Table 2 (FCQ) and Table 3 (MPL). For radii < 0.7" from the center 
the data were sampled at full resolution(~ 0.05") while at larger radii they were binned 
spatially. The sampling at small radii is fine enough that the data points are somewhat 
correlated; this was done to ensure that no information concerning the steep radial gradients 
of the profiles was lost. 

3.2. Results for M32 

Figure |] presents LOSVDs computed via the MPL algorithm at positions separated 
by about 0.1" along the slit. One expects these broadening functions to obey 
N(V;R) = N(—V; —R), since for a point-symmetric galaxy, the velocity distributions 
should reverse after passing from one side of the galaxy to the other. The LOSVDs 
of Fig. |] show approximately the expected symmetry. The right-hand column of 
Figure ^ plots mean broadening functions averaged over the two sides of the galaxy, 
N(V) = \\N{V, R) + N(-V, -R)}; the central LOSVD has been symmetrized about V = 0. 

These broadening functions show clear and consistent deviations from Gaussian form, 
in two respects. First, the central LOSVD exhibits strong super-Gaussian "wings" at high 
velocities. These wings are possibly present also in some of the off-center LOSVDs although 
with lower amplitude. Second, the off-center LOSVDs are asymmetric, with tails extending 
toward velocities opposite in sign to the mean velocity at each radius. These asymmetric 
tails are similar to those exhibited by a rotating system superimposed on a slowly-rotating 
bulge. 

Systematic problems in the spectral deconvolution, e.g. template mismatch or incorrect 
continuum subtraction, can easily produce features like the wings and tails seen in the 
broadening functions of Figure [|. Such errors in most cases would be expected to produce 
features located at the same velocity on both sides of the galaxy (e.g. [Bender, Saglia fc| 
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Gerhard 1994) and are therefore an unlikely explanation for the asymmetric tails seen in 
the off-center LOSVDs. The strong wings seen in the central LOSVD might more plausibly 
be attributed to systematic errors. However we found that the wings in the central LOSVD 
were robust; they appeared in both the MPL and FCQ estimates of N(V) (though less 
clearly in the latter - see Figure and were relatively unaffected by changes in the assumed 
continuum level or slope. We carried out MPL deconvolutions where the fit to the galaxy 
spectrum was restricted to the region around only one, or two, of the three calcium-triplet 
lines; these LOSVDs also exhibited strong wings. We also tried using one of the other 
available STIS stellar templates; again the wings were only slightly affected. (The adopted 
template, Figure 1, produced the best overall fit to the galaxy spectrum.) Finally, we show 
in §4.2 that the wings are consistent with those predicted by stars in the gravitational field 
of a super massive black hole. 

Figure |5] shows the Gauss-Hermite parameter Vq, a measure of the stellar rotation, 
in the inner arc second of M32. Also plotted are /13, the lowest, odd GH moment of 
the LOSVD, and the "corrected" rotation velocity, V 0iC = V + v^3o"o^3- Vo tC is a closer 
approximation than Vq to the true mean line-of-sight velocity V ( |van der Marel fc Franx 



1993|) . The corrected rotation velocity is lower in absolute magnitude than |Vq| due to the 



asymmetric wings of the LOSVD noted above. 

The STIS rotation curve is consistent with earlier ground-based measurements (Figure 
0) at radii > 1" but with a larger peak value, ~ 60 km s -1 . Furthermore the rotation 
curve remains flat or slightly rising into smaller radii than seen heretofore, before falling 
at R < 0.25" due to the blending of light from the two sides of the galaxy. There is a 
suggestion of an east-west asymmetry in the rotation curve though the effect is probably 
not significant. 

The hs profile is approximately antisymmetric about the center of M32, as expected 
in a relaxed galaxy. \h^\ reaches a maximum value of ~ 0.05 at \R\ ~ 0.3" and appears to 
gradually decline at larger radii. This behavior is similar to that predicted in axisymmetric 
models (e.g. Figure 8 of |Dehnen 1995|) where /13 remains essentially constant at radii outside 
the seeing disk. 

The Gauss-Hermite parameter <7o is shown in Figure ^. Also plotted is the lowest, 
even moment of the LOSVD, and the "corrected" velocity dispersion, ct 0jC = cr (l + \f^h^)\ 
o"o jC is a closer approximation than oq to the true rms velocity a. The velocity dispersion 
rises suddenly inside of ~ 0.3", approximately the same radius at which the rotation curve 
begins to fall. This coincidence suggests that at least part of the rise in a is due to 
averaging of the rotation velocity over the two sides of the galaxy near the center, which 
has the effect of converting a rotation into an apparent dispersion QTonry 1987|) . 



-14- 



The corrected velocity dispersion o"o >c rises well above o"o near the center due to the 
strong non-Gaussian wings of the LOSVD. The central value of o"o jC is ~ 175 km s _1 ; this 
should probably be interpreted as a lower limit since /14 is only sensitive to the inner parts 
of the wings. (We argue below, based on model fits, that the rms velocity in the central 
resolution element may be as high as ~ 200 km s _1 .) The ground-based data (Figure [?]) are 
consistent with the STIS dispersions at radii > 1" but fail to resolve the continued rise in 
cr inside of ~ 0.5". 



Dynamical models (e.g. Pehnen 1995| ; |Qian et al. 1995| ) predict h^{R) profiles similar 



to that in Figure [!] when observed with HST resolution: a central maximum; a rapid drop, 
to small or negative values, at R ~ 0.1"; and a nearly constant value at larger radii. The 
predicted drop at ~ 0.1" is due to blending of the light from the two sides of the galaxy, 
which broadens the low-velocity part of the LOSVD and lowers the observed h^. The 
predicted central value of h± depends strongly on the black hole mass and on the PSF; our 
value, hi ~ 0.14, is larger than in the two studies just cited, but these studies were based 
on rather low assumed black hole masses, Mh = 1 — 2 x 10 6 Ai Q . The true black hole mass 
is probably greater (van der Marel et al. 199S ) . 



The behavior of at larger radii is surprising. Previous observational studies (e.g. van 



der Marel et al. 1994a] ; [Bender, Kormendy fc Dehnen 19961) have returned smaller estimates 



for /i 4 in M32. However we believe that these earlier results are not inconsistent with ours 
giving the difficulties involved with estimating this parameter. The van der Marel (1994a) 
study was based on WHT observations with a much lower spatial resolution than the STIS 
data. At radii < 1", the value of recovered by those authors was strongly affected by the 
PSF blending of the rotation curve discussed above, yielding negative values in the central 
aperture. Outside of ~ 2", van der Marel et al. found h^ to increase sharply to ~ 0.03 on 
both sides of the galaxy (their Figure 12). Van der Marel (private communication) notes 
that the values of /14 derived from the WHT data depended sensitively on the choice of 
template spectrum and on the algorithm for continuum subtraction. Using a single, best-fit 
template, was found to lie between ~ 0.3 and ~ 0.5 throughout the inner 2"; the lower 
values of in the published paper were derived using a spectral deconvolution routine 
that constructs an "optimal" template by linear superposition of a set of stellar spectra. 
In another ground-based study, Bender, Kormendy & Dehnen (1996) applied the FCQ 
algorithm to CFHT data of higher spatial resolution and found « 0.05 inside of 0.2", 
gradually falling to ~ at ~ 1.0". However the spectral resolution in this study was only 
80 km s _1 and the derivation of h 3 and correspondingly difficult; as noted above, we also 
found hi m from the STIS data using the FCQ algorithm and argued that these values 
were significantly negatively biased. 
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Although we believe that all of these studies are consistent with our conclusion that 
/i 4 is significantly positive throughout the nucleus of M32, we are less willing to strongly 
endorse the precise values shown in Figure due to the sensitive dependence of this 
parameter on the details of the spectral deconvolution algorithm, continuum subtraction, 
smoothing level, etc. In Paper II we will construct dynamical models based on a range of 
assumed profiles in order to test the dependence of the inferred black hole mass on this 
parameter. 

We may also compare our results to the van der Marel et al. (1997, 1998) HST/FOS 
measurements of Vq and o"o (Figure §). The FOS measurements were taken through square 
apertures as small as ~ 0.1" on a side, hence their spatial resolution is comparable to 
that of the STIS data. However the FOS is a low spectral resolution instrument and not 
well suited to objects like M32 with a relatively low velocity dispersion; furthermore there 
are difficulties in positioning the FOS and these were probably the cause of the large 
point-to-point variations seen by van der Marel et al. One advantage of STIS over FOS is 
the continuous spatial sampling which avoids potential errors in aperture placement. We 
find a hint in the STIS data of the asymmetry seen in the FOS <Tq(R) profile (a more rapid 
falloff on the west side). The central FOS value of <To seems significantly bigger than found 
here, and the FOS rotation velocities are systematically lower. 
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4. Analysis 

The STIS data show what appear to be clear signatures of the gravitational influence 
of a massive compact object on the stellar velocity distribution within the central parsec 
of M32. Here we use simple dynamical models to address the question of whether these 
observed features are consistent with the presence of a supermassive black hole. Our aim 
is not to derive the best possible estimate of the black hole mass or the stellar velocity 
distribution - those are the goals of Paper II - but rather to address two, more basic issues 
concerning the interpretation of the data. 

1. The velocity dispersion profile (Figs. |6], 0) exhibits a sudden upturn at a distance 
of ~ 0.3" from the center, presumably due in part to the gravitational force from a 
massive compact object. At roughly the same radius, the rotation curve falls (Figs. |], |7|), 
presumably due to blending of light from opposite sides of the galaxy which are rotating in 
opposite directions. The blending would also be expected to contribute to the rise in the 
observed velocity dispersion ( |Tonry 1987|) , consistent with the fact that the upturn in a 
and the drop in V occur at roughly the same radius. We would like to estimate the degree 
to which the velocity dispersion spike is a product of this blending, and the degree to which 
it is due to a real upturn in the stellar random velocities. In other words: Do the STIS data 
resolve the black hole's sphere of influence? 

2. The central LOSVD in M32 exhibits strong, super-Gaussian wings (Fig. |2|, ^). 
Such wings are a generic prediction of the black hole model ( Bahcall fc Wolf 1976 



van aer 



Marel 19"94] ; Pehnen 1995|) ; they result from stars on high- velocity orbits within the black 



hole's sphere of influence. However systematic errors in the spectral deconvolution can also 
produce spurious features in the LOSVD's, particularly at large velocities where the form 
of the broadening function is only weakly constrained by the spectra. We would like to 
verify that the inferred, non-Gaussian shape of the central LOSVD is consistent with that 
expected from the black hole model. 

We note that M32 is the only galaxy so far to exhibit either a resolved central spike in 
the stellar velocity dispersions, or strong non-Gaussian wings in the nuclear LOSVD. Either 
feature, if observed with sufficient spatial resolution and S/N, could independently place 
strong constraints on the mass of a central black hole. We will in fact generate estimates of 
Mh from our analyses of both features but we stress that the best estimates of Mh can only 
come from more complete modelling based on the entire kinematical data set. 
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4.1. The Velocity Dispersion Spike 



One expects to see a rise in the stellar velocities in a hot stellar system at a distance 
~ r/j from the black hole, where 

GM h 

r h = —3-, (2) 

* 

the black hole's "radius of influence" ( [Peebles 1972|) . Here is the black hole mass and 



a 2 is the stellar velocity dispersion at r > r^. Setting = 10 6 .M© and a* = 60 km 
s _1 (Figures |7]) gives ~ 1.2 pc ~ 0.3". This value seems comfortably larger than 
the HST/STIS resolution but it is based on an assumed value of Mh and furthermore it 
refers to the true radius whereas we observe the galaxy in projection, which tends to hide 
otherwise sharp features. Modelling the spike therefore requires us to predict the true, 
projected velocity field of the galaxy in two dimensions on the plane of the sky, including 
both random and rotational velocities, and then to convolve it with the STIS PSF. 

We begin by constructing solutions to the stellar hydro dynamical equations. We 
assume the simplest possible axisymmetric model for the stars, in which the velocity 
dispersions are isotropic in the meridional plane (w, z), i.e. o ro (ca, z) = <J z (zu, z) = <j(zu, z); 
the model is flattened by an inequality between a 2 and the mean square azimuthal velocity 
v^. The second moments of the internal velocity distribution are given by 

roc d$ 5$ d(ua 2 ) 

vo = / u—dz, uvi = va + vzu- V w — (3) 

Jz oz v ozu ozu 

where v is the stellar number density, is the azimuthal velocity and $ is the combined 

gravitational potential from the stars and the central black hole, Q(zu, z) = $*(ro, z)—GMh/r 

( Merritt 1999| ). 



We evaluated these expressions assuming a stellar density 

u(zu, z) = u (m/b) a [l + (m/b) 2 ] [l + (m/c) 2 ] 7 , m 2 = zu 2 + (z/q) 2 , (4) 

a parametrized form proposed by van der Marel et al. (1998); those authors 
found a good match between their model (observed edge-on) and M32 with the 
parameters a = -1.435, (3 = -0.423,7 = -1.298,6 = 0.55", c = 102.0", q = 0.73 and 
uo = 0.463 x lO 5 L 0i ypc -3 . The luminosity density uq is converted into a mass density po by 
the factor (M/L)(A4q/Lq), with M/L the mass-to- light ratio of the stars in solar units. 

Given a(zu,z) and v 2 (zu,z) as obtained from equations (3) and (4), the projected, 
line-of-sight, mean square velocity V 2 is obtained by a density- weighted integration through 
the galaxy, assumed here to be edge-on: 

R 2 



POO 

E(R,Z)V 2 (R,Z) = 2 / u(zu,z 

Jr 




a 2 (zu,z) + ^vl(zu,z) 
zu z v 



zudzu 
Vzu* - R* (5) 
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( |Fillmore 1986|) , where (R, Z) are coordinates on the plane of the sky and R is measured 
parallel to the long axis of the galaxy's figure; S(i2, Z) is the stellar surface density. 

Values of V 2 were computed on a rectangular grid of 180 x 25 locations with separations 
of 0.015" in R and 0.02" in z. These values were then convolved with the STIS PSF 
and averaged over the pixel area and the aperture after weighting by the model surface 
brightness. The STIS PSF at 8500A has a FWHM of ~ 0.115". The PSF is also slightly 
asymmetric ( Power et al. 20001) . We ignored this slight asymmetry when carrying out the 
convolutions with our models. 

The second velocity moments of models constructed in this way are uniquely determined 
by the two parameters (M^, M/L) that specify the potential. Figure |9| shows the goodness 
of fit of the models to the observed, mean square velocities; only data points within the 
inner 1.0" were used in evaluating x 2 ■ As estimates of V 2 , we took Vq c + a 2 c , where V 0tC 
and <To, c are the Gauss-Hermite parameters corrected by h 3 and respectively (Figures 
|, |). The best-fit model has M h « 3.2 x 10 6 .M Q and M/L « 3.3 with X 2 = 0.64; a x 2 of 
unity includes models with Mh as small as 2.2 x 10 6 .Mo and as large as 4.3 x 10 6 Ai & . 

The degree of net rotation in our models may be adjusted by partitioning the azimuthal 
motions between streaming, v$, and dispersion, a^. We followed the standard practice 
( [Satoh 1980| ) of making <j\ a weighted average of a 2 and v\, i.e. 



a 2 , = k 2 a 2 + (1 - k 2 )v 2 . (6) 

The parameter k (assumed independent of position) may be varied between zero 
(corresponding to no streaming motions) and a maximum value, of order unity, at which a 2 , 
is forced below zero at some point in the meridional plane; k — 1 yields an "isotropic oblate 
rotator." 

When we add A; as a free parameter and require the models to fit the observed rotation 
and velocity dispersion profiles separately, the best-fit values of (Mh, M/L) were nearly 
unchanged but x 2 increased to 3.7 - since the model is now being asked to fit twice as many 
data points with only one extra parameter. The best-fit value of k was found always to be 
close to 1.2, implying slightly smaller a§ (i.e. greater rotation) than in an isotropic oblate 
rotator. Figure [1(] compares the data with the predicted profiles for Mh = 3.0 x 1O 6 A^0, 
close to the best-fit value, and for M h = 2.0 and 4.0 x 1O 6 A^ . 



We draw the following conclusions from these comparisons. 

1. The lowest-order moments of the line-of-sight velocity distribution in M32 are 
reasonably well fit near the center by our simple axisymmetric model, with a black hole 
mass M h ^3x 1O 6 .M . The rotation curve is best fit by a smaller mass (~ 1 — 2 x 10 6 .M ©) 
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and the velocity dispersions by a larger mass (~ 3 - 4 x 10 6 A4 & ); if we require the models 
only to fit the mean square velocities, the fit is essentially perfect within the inner arc 
second. 

2. The STIS observations probably come close to resolving the central upturn in the 
stellar velocity dispersions, which is predicted to occur at a projected radius of ~ 0.1" for 
M fc ~2x 10 6 M Q and ~ 0.2" for M h ~ 4 x 1O 6 .M . 

3. Smearing of the stellar rotation field probably accounts for only a small part of the 
observed upturn in the dispersions. 

We note that Dehnen (1995) was able to improve the fit of his axisymmetric models to 
the M32 data then available by varying the ratio of rotational to non-ordered azimuthal 
motions, which in our models would correspond to varying k with position. 

If our simple model for the internal dynamics of M32 is approximately correct, Figures 
[| and imply that the STIS data can place strong constraints on the mass of the 
central black hole. However such a conclusion must await the results of the more complete 
modelling of Paper II. That study may yield tighter constraints on Mh, due to the use 
of the full kinematical data set; or weaker constraints, due to the increased flexibility of 
general, anisotropic models. 



4.2. The Central Broadening Function 

The LOSVD near the projected center of a galaxy containing a black hole is expected 
to be very non-Gaussian due to high velocity stars orbiting near the central mass (|Bahcall~SE| 
Wolf 1976|) . If the stellar density follows a power law into the center, p oc r -7 , the LOSVD 



in an aperture containing the black hole has power-law wings, N(V) oc V 2 " 1 7 , V — > oo ( Ivan 



der Marel 1994| ). The amplitude of these wings depends on the ratio of black hole mass to 



stellar mass within the aperture, and on the slope 7 of the stellar density profile ( |Dehnen 



1995), among other factors. The wings are most prominent in a galaxy, like M32, for which 
the stellar cusp is steep, 7 ~ 2, since a large fraction of the light near the projected center 
comes from the region near the black hole. 

Here we ask whether the strong wings seen in the central M32 broadening function 
(Figures |2], |j) are consistent with the black hole model. To answer this question we 
must compute a stellar distribution function / and integrate it over the two velocity 
components in the plane of the sky ( |Merritt 1987] ) . The result must then be smeared by 
the instrumental resolution. We once again restrict ourselves to the simplest model which 
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permits a meaningful comparison with the data; in this case, a spherical galaxy with an 
isotropic distribution function, / = f(E). Our model will ignore rotation, which acts to 
broaden the central LOSVD but leaves it symmetric. In any case we do not know what the 
contribution of rotation to the stellar velocity distribution is very near to the black hole. 
For the stellar density profile we assume the spherical version of equation (4); f(E) then 
follows from the standard formula ( [Eddington 1916|) and the projected velocity distribution 
is also straightforward to compute flMerritt 1993| ). We chose to fix M/L for the stars at 
3.5(.M©/ L Q ) as the black hole mass was varied; this M/L reproduces approximately the 
observed mean square velocities for r ^> r^. 

Figure [XT] compares the central LOSVD in M32 to the broadening functions predicted by 
the spherical model, for black hole masses M h = (2.5, 5.0, 10.) x 10 6 M Q . The M32 LOSVD 
has been symmetrized about V = 0. The fit is quite reasonable for Mh = 5.0 x 10 6 .M©, 
although the high-velocity wings are better fit by still larger masses. The model LOSVD 
for Mh = 5.0 x 10 6 A4 Q has the parameters: 

a = 178 km/s a = 133 km/s h 4 = 0.085. (7) 

The first two numbers are essentially identical to the values inferred from the M32 data 
(Figure |6|); /14 is lower than, but consistent with, the M32 estimate (0.14 ±0.03). 

We conclude that the central LOSVD in M32 is consistent with that expected for a 
stellar nucleus containing a massive compact object, with a mass comparable to that found 
in the fit to the axisymmetric models. 
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5. Summary 

We used HST and STIS to obtain stellar absorption line spectra near the center of 
M32 in a wavelength region centered on the Calcium triplet. The spectra were analyzed 
using two independent spectral deconvolution routines; these gave fully consistent results 
except in the case of the Gauss-Hermite /14 parameter, but we argued that the differences 
could be reconciled after taking into account the different biases of the two algorithms. 
The stellar rotation velocities in M32 are slightly higher than observed from the ground 
and remain constant into ~ 0.25" from the center. The velocity dispersions exhibit a clear 
spike beginning at approximately the same radius. These two kinematical profiles are 
consistent with those predicted by simple axisymmetric models containing central black 
holes with masses in the range 2 — 5 x 10 6 .M©. The stellar LOSVDs show significant 
deviations from Gaussian form as measured by the Gauss-Hermite parameters h 3 and /14. 
The central LOSVD is particularly non-Gaussian, exhibiting strong, high-velocity wings. 
We showed that the amplitude of these wings is consistent with that predicted by simple 
models containing black holes with masses of order 3 x 10 6 .M©. 

Detailed dynamical modelling of M32 based on these data and estimates of the black 
hole mass will be presented in Paper 11. 

We thank W. Dehnen and R. van der Marel for helpful discussions. This work was 
supported by NASA grants NAG 5-3158 and NAG 5-6037, by NSF grant AST 96-17088, 
and by STIS GTO funding. Data presented here were based on observations with the 
NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in Astronomy, Inc. 
(AURA), under NASA contract NAS5-26555. 
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Appendix A 
Gauss-Hermite Moments 

The two spectral deconvolution algorithms described above yield nonparametric 
estimates N(V) of the stellar LOSVD. Here we describe the methods used by the two 
algorithms to derive the GH moments from N(V). 

Let N(X, Y; V) be the distribution of line-of-sight stellar velocities in the aperture 
centered at (X, Y). Define the GH moments of N as 

/oo 
N(X,Y;V)g(w)Hi(w)dV, (Al) 
-oo 

where Hi are the Hermite polynomials (as defined by |Gerhard 1993|) and the weight function 

g(w) = —^e~ w2 / 2 , w = (V- V )/a (A2) 

V27T70 

has three free parameters (70, Vq, <To). Following van der Marel Sz Franx (1993), we choose 
these parameters at every point (X, Y) such that 

h (X,Y) = l, h 1 (X,Y)=h 2 (X,Y)=0. (A3) 

These definitions impose the following implicit conditions on (70, Vo, cr ): 

70 = V2a l°° N{V)e- w2/2 dw, (A4a) 

/oo „ 
N(V)e~ w 12 w dw } (A4b) 
-00 

roc , 

= / N(V)e~ w /2 {2w 2 - 1) dw. (A4c) 



The relations (A4) define a nonlinear minimization problem with solutions (70, Vo, cr ) given 
N(V). 

The MPL algorithm ( [Merritt 1997|) derives the three parameters in just this way, 
using the NAG routine E04FDF to minimize the sum (ho — l) 2 + h\ + h\ as a function 
of (70, Vo, (To). The higher-order GH moments are then derived using equation (Al), by 
numerical integration over N(V). 

Most spectral deconvolution algorithms of which we are aware derive the parameters 
(70, Vq, (To) in a different way. The LOSVD is compared to the trial function 



v2vro"o 



J max 



l+5>^» 

J'=3 



(A5) 
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where j max is the index of the highest GH moment fitted to N(V); typically j max = 4. The 
integrated square deviation between N(V) and A/"(V) is then minimized by varying the 
(jmax + 1) free parameters (■j ,V ,a ,h 3 ,h 4: ,..., hj max ). This is the technique used by the 
FCQ algorithm. 



A theorem QMyller-Lebedeff 1908|) guarantees the equivalence of the two approaches if 



jmax = oo in equation (A5) ([van der Marel fc Franx 1993|) . However if j max ^ oo, and if 
the input N(V) can not be precisely represented by a finite GH series with j < j max , the 
results given by the two algorithms will differ. For instance, in attempting to represent an 
N(V) having he ^ using j max = 4, the FCQ algorithm will adjust o"o and h± to incorrect 
values in order to better fit the high-velocity wings of the profile with the limited number 
of terms allowed to it. This is illustrated in Figure Al, which shows the values of o"o and 
hi generated by the second algorithm, <7 and h\, compared to the true values for an input 
N(V) with hi = 0.15 and nonzero h^. 

N(V) = ^=e- y2 l 2 {1 + 0.15F 4 (V) + h 6 H e (V)} . (A6) 
v 2tt 



For \he\ > 0.1, the errors in <7o and h± as derived from the second algorithm are > 15% and 
> 20% respectively. 
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Appendix B 

Performance Evaluation of the FCQ and MPL Algorithms 

Here we compare the performance of the FCQ and MPL algorithms given simulated 
data. Our primary goal is to understand the source of the systematic offset of h A values as 
derived from the M32 spectra by the two algorithms (§3.1). Two independent sets of tests 
were carried out, the first by R. Bender and the second by D. Merritt. All tests were based 
on synthesized galaxy spectra generated from the STIS template spectrum (Figure [I]) by 
convolving it with an assumed N(V) and adding noise. 

The first set of tests addressed the accuracy of FCQ estimates when the galaxy velocity 
dispersion is low. It is well known that the accuracy of FCQ begins to fall off when the 
galaxy velocity dispersion becomes comparable to the dispersion of the template spectrum 
(e.g. Pender, Paquet fc Nieto 1991|) . Figure Bl shows values of <Jo recovered by FCQ 
given a Gaussian-broadened template spectrum and thirty random noise realizations. The 
default value (W ~ 1) of the smoothing parameter was used. There is a positive bias in 
the estimated values beginning at cr ~ 100 km s _1 ; the bias increases with decreasing 
do becoming significant for <To ~ 50 km s^ 1 . The bias is only weakly dependent on S/N. 
This bias in the estimation of <jq is unlikely to be important for the nucleus of M32 where 
cT Z 100km s -1 . 

Figure B2 shows the performance of FCQ at recovering /i 4 . The template spectrum 
was broadened using an N(V) of the form 

Ni(V) = -±=e- y2 l 2 °» {1 + h A H A {V/a Q )} (Bl) 
V In 

with various values of o"o and h A . Figure B2 reveals significant biases in h A for <jq < 100 km 
s _1 , even when S/N is as great as 100. When <jq ~ 50 km s -1 and S/N m 30, characteristic 
of M32 at ~ 1", the bias in h A is ~ —0.1 for an input h A of ~ 0.1. 

The second set of tests compared the performance of the FCQ and MPL algorithms on 
galaxy spectra generated from the broadening function 

N *W = — TXTTTT^' (B2) 
no 1 + [y joy 

a Lorentzian function, with a = 100 km s _1 . This LOSVD is qualitatively similar to what 
is expected in a black-hole cusp, with iV ~ V~ 2 high-velocity wings. The non-trivial GH 
parameters are 

7o = 0.76986 a = 108.07 km s _1 h± = 0.14546 h 6 = 0.01850. (B3) 
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Figure B3 shows mean estimates of N(V) obtained using the two algorithms for 100 random 
realizations of the noise. The smoothing parameter in both algorithms was adjusted to 
minimize the mean square error of N(V) (as defined below) for each value of S/N. There is 
a greater bias in the FCQ estimates, as well as a persistent "ringing" at high velocities. 

Figure B4 plots the mean integrated square error (MISE) and integrated square 
bias (ISB) of the recovered broadening functions as functions of S/N; in the case of the 
FCQ algorithm, the integrated errors are shown both for the optimal choice of smoothing 
parameter W op t that minimizes the MISE, as well as for the value chosen by the algorithm 
(~ 1.3). The MISE of an estimate f(x) is defined as 

MISE [f(x)] = E J [f(x) - f(x)} 2 dx (B4a) 

= J { E f(x)-f(x)} 2 dx + J (E{f 2 {x)}-E{f{x)f^dx (B4b) 

= ISB [f(x)] + IV [f(x)] , (B4c) 

the sum of the integrated square bias ISB and the integrated variance IV ( |Silverman 1986| ); 
here E denotes the expectation value, i.e. the average over many random realizations of 
the noise. The MISE and ISB displayed in Fig. B4 were divided by the normalizing factor 
J[N 2 (V)} 2 dV. 

The MISE of the MPL estimates falls roughly as a power law, MISE[iV(V)] ~ (S/N)' 1 - 3 , 
close to the asymptotic (S/N)' 1 of parametric estimators. Approximately 1/2 of the total 
square error comes from the bias and 1/2 from the variance. In the case of the FCQ 
algorithm, the MISE behaves in a more complicated way with S/N, at first falling with 
S/N then appearing to level off for S/N > 50. This levelling off is a consequence of the 
low- velocity-dispersion bias of FCQ discussed above. For the FCQ estimates, the bulk of 
the MISE is due to the variance; adjusting the smoothing parameter primarily affects the 
bias and has little effect on the MISE. For S/N w 20, the mean square error of the optimal 
FCQ estimate is a factor ~ 3 greater than that of the MPL estimates. 

The bias in N(V) is in the direction of wider and more steeply truncated functions, 
particularly in the case of the FCQ estimates (Fig. B3). This bias in N(V) is consistent 
with the negative bias found above in estimates of h±. Figure B5 compares the ability of 
the two algorithms to recover from the Lorentzian A^(V). Plotted there are the mean 
square error (MSE) and bias in estimates of from 100 random noise realizations; the 
MSE is defined, for any estimated parameter P, as 

MSE(P) = E {P - P} 2 , (B5) 
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which can also be decomposed into contributions from the squared bias SB and the variance 
V: 

MSE(P) = (EP -Pf +{e {P 2 } - E {P} 2 ^j (B6a) 

= SB(P) +V(P). (B6b) 

The MSE of estimates obtained with the MPL algorithm again varies roughly as a power 
law, MSE(h 4 ) ~ (S/N)^ 1 - 5 . The bias in the MPL estimates is always negative, i.e. in the 
direction of more Gaussian JV(V)'s; for S/N ~ 20, this bias is a modest ~ —0.03, dropping 
to < -0.01 for S/N = 100. 

The FCQ estimates of h 4 show a considerably greater error, both in the bias and 
the variance. Two sets of FCQ estimates were made: first using the default value of the 
smoothing parameter returned by the code, W ~ 1.3; and second using the optimum value 
Wopt that minimized the MSE of the h 4 estimates at each S/N. For the default value of W, 
the algorithm returns mean estimates of h 4 that lie in the range 0.03 — 0.05 for all values 
of S/N > 10, an average error of ~ 70%. However the optimum smoothing parameter for 
the recovery of h 4 was found to vary strongly with S/N, from W op t ~ 0.5 for S/N = 5 to 
W^t ~ 2 for S/N = 100 (Fig. B5c). Nevertheless a substantial bias remains when W op t is 
used, of order ~ —0.05 even for S/N = 50 — 100. These biases are larger than found above 
using a more Gaussian N(V) with smaller h 4 and suggest that FCQ estimates of h 4 may be 
substantially biased even for cr as large as ~ 100 km s^ 1 when the true N(V) is sufficiently 
non-Gaussian. 
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Table 1. Observational Setup 



Gain (e"/ADU) 


1.0 


Wavelength coverage 


8275A - 8847A 


Reciprocal dispersion (A pixel -1 ) 


0.56 


Aperture 


52" x 0.1" 


Comparison line FWHM (pixels) 


2.0 


R = A/AA 


7644 


Instrumental dispersion (07) (km s -1 ) 


17.1 


Spatial scale (" pixel -1 ) 


0.05071 
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Table 2. M32 Kinematics as Derived via FCQ 



K 


T T 


AT/ 

Al/ 


O"0 




h 3 


A U 


h,4 


A U 


T T 

V ,c 


AT/ 


O"0,c 




I 1 ) 


(2) 




(3) 




(4) 




(5) 




(6) 




/7\ 
( 7 ) 




4.497 


25.2 


32.2 


67.2 


19.0 


0.110 


0.436 


-0.289 


0.436 


38.0 


60.2 


19.6 


75.5 


4.000 


49.7 


31.5 


97.5 


32.2 


0.053 


0.294 


-0.055 


0.294 


58.7 


58.9 


84.3 


77.3 


3.465 


55.5 


23.6 


60.7 


13.9 


-0.063 


0.353 


-0.230 


0.353 


48.8 


44.0 


26.5 


54.9 


2.983 


54.4 


20.6 


82.2 


18.9 


0.017 


0.227 


-0.090 


0.227 


56.8 


38.3 


64.1 


49.6 


2.546 


35.0 


18.8 


60.1 


14.0 


0.040 


0.284 


-0.148 


0.284 


39.1 


35.0 


38.3 


44.4 


2.173 


24.1 


25.1 


103.3 


18.7 


-0.052 


0.221 


-0.148 


0.221 


14.8 


46.9 


65.9 


59.4 


1.842 


35.9 


17.2 


64.0 


13.0 


-0.058 


0.245 


-0.144 


0.245 


29.5 


32.2 


41.4 


40.8 


1.562 


61.5 


17.2 


88.8 


15.6 


-0.074 


0.176 


-0.093 


0.176 


50.1 


32.2 


68.6 


41.5 


1.336 


42.8 


13.7 


76.9 


13.5 


-0.013 


0.162 


-0.067 


0.162 


41.0 


25.6 


64.3 


33.4 


1 1 X.Q 




10.41 


81 P. 
ol.O 


1 3 Q 

lo.y 


n nm 

U.U1U 


n 1 71 

U. 1 ( 1 


-u.uyo 


n 1 71 

U.l 1 1 


OO.O 


98 P. 




^7 n 


1.034 


53.6 


15.3 


82.9 


12.6 


0.059 


0.168 


-0.122 


0.168 


62.1 


28.6 


58.1 


36.6 


0.932 


52.8 


13.5 


74.4 


11.2 


-0.024 


0.165 


-0.120 


0.165 


49.7 


25.2 


52.5 


32.3 


0.831 


65.7 


12.9 


72.2 


10.4 


-0.078 


0.163 


-0.129 


0.163 


55.9 


24.2 


49.4 


30.8 


0.729 


69.6 


11.2 


77.1 


11.3 


-0.001 


0.132 


-0.057 


0.132 


69.4 


20.9 


66.3 


27.4 


0.654 


62.4 


12.6 


72.5 


11.0 


0.018 


0.158 


-0.105 


0.158 


64.7 


23.5 


53.9 


30.3 


0.603 


57.5 


12.6 


79.0 


8.1 


-0.055 


0.144 


-0.180 


0.144 


50.0 


23.4 


44.2 


29.2 


0.553 


56.8 


13.0 


75.3 


11.1 


-0.014 


0.157 


-0.112 


0.157 


55.0 


24.3 


54.7 


31.2 


0.502 


58.8 


10.7 


72.4 


10.1 


0.019 


0.135 


-0.083 


0.135 


61.1 


20.1 


57.7 


26.1 


0.451 


58.8 


11.6 


81.1 


11.2 


0.000 


0.130 


-0.074 


0.130 


58.8 


21.6 


66.4 


28.2 


0.401 


65.2 


10.4 


76.0 


10.2 


-0.089 


0.124 


-0.065 


0.124 


53.5 


19.4 


63.9 


25.3 


0.350 


58.8 


10.0 


86.0 


10.2 


-0.067 


0.106 


-0.055 


0.106 


48.8 


18.7 


74.4 


24.6 


0.299 


52.9 


9.6 


89.3 


10.3 


-0.072 


0.098 


-0.037 


0.098 


41.8 


18.0 


81.2 


23.8 


0.248 


60.6 


8.4 


77.4 


8.9 


-0.041 


0.099 


-0.044 


0.099 


55.1 


15.7 


69.1 


20.8 



Table 2 — Continued 



K 


T T 

Vo 


AT/ 


O"0 


A 

A(T 




A U 


h,4 


A U 


Vq,c 


AT/ 


O"0,c 




I 1 ) 


(2) 




(3) 




( 4 ) 




(5) 




(6) 




( 7 ) 




0.198 


53.4 


8.1 


90.4 


8.3 


-0.049 


0.082 


-0.057 


0.082 


45.7 


15.2 


77.8 


20.0 


0.147 


40.8 


9.2 


99.4 


8.9 


0.038 


0.084 


-0.072 


0.084 


47.3 


17.2 


81.8 


22.4 


0.096 


38.9 


10.2 


112.2 


11.2 


0.012 


0.083 


-0.030 


0.083 


41.2 


19.1 


104.0 


25.4 


0.046 


36.3 


8.7 


123.0 


10.8 


-0.009 


0.064 


0.019 


0.064 


34.4 


16.2 


128.7 


22.1 


-0.005 


-2.9 


9.0 


136.6 


12.6 


-0.013 


0.060 


0.075 


0.060 


-6.0 


16.8 


161.7 


23.8 


-0.056 


-24.6 


9.6 


124.2 


10.1 


0.032 


0.071 


-0.045 


0.071 


-17.7 


18.1 


110.5 


23.9 


-0.106 


-34.4 


9.2 


104.2 


9.8 


0.070 


0.081 


-0.041 


0.081 


-21.7 


17.3 


93.8 


22.9 


-0.157 


-40.9 


9.9 


98.3 


11.5 


0.087 


0.092 


-0.007 


0.092 


-26.1 


18.6 


96.6 


25.0 


-U.ZUo 


-OO. 1 


Q Q 

y.y 


on 8 
yu.o 


i n a 


n nQ7 
U.Uo ( 


n 1 nn 

U. 1UU 


-U.U40 


n 1 nn 

U. 1UU 


AT Q 


lo.O 


sn 8 
oU.o 


OA P. 


-0.259 


-64.0 


9.9 


87.3 


10.6 


0.013 


0.103 


-0.036 


0.103 


-62.0 


18.4 


79.6 


24.5 


-0.309 


-62.6 


9.0 


76.4 


8.5 


0.052 


0.107 


-0.079 


0.107 


-55.7 


16.8 


61.6 


21.8 


-0.360 


-63.4 


10.2 


80.8 


10.8 


0.045 


0.114 


-0.042 


0.114 


-57.1 


18.9 


72.5 


25.0 


-0.411 


-50.9 


12.0 


90.4 


13.8 


0.038 


0.120 


-0.010 


0.120 


-44.9 


22.3 


88.2 


30.0 


-0.461 


-57.6 


11.6 


70.5 


11.0 


0.016 


0.150 


-0.078 


0.150 


-55.7 


21.7 


57.0 


28.2 


-0.512 


-60.1 


10.8 


63.5 


9.0 


0.033 


0.155 


-0.118 


0.155 


-56.5 


20.2 


45.1 


25.9 


-0.563 


-55.0 


13.5 


83.4 


12.1 


0.092 


0.147 


-0.096 


0.147 


-41.7 


25.2 


63.8 


32.5 


-0.614 


-51.5 


15.5 


93.2 


15.2 


0.047 


0.151 


-0.067 


0.151 


-43.9 


28.9 


77.9 


37.8 


-0.664 


-51.4 


13.4 


66.4 


11.6 


0.053 


0.183 


-0.106 


0.183 


-45.3 


25.0 


49.2 


32.1 


-0.715 


-58.4 


13.2 


78.0 


10.7 


0.173 


0.154 


-0.126 


0.154 


-35.0 


24.8 


53.9 


31.5 


-0.766 


-46.2 


18.0 


87.3 


17.3 


0.189 


0.188 


-0.076 


0.188 


-17.6 


34.1 


71.0 


43.9 


-0.841 


-31.8 


12.2 


77.9 


10.1 


-0.079 


0.142 


-0.119 


0.142 


-42.4 


22.7 


55.2 


29.1 


-0.943 


-51.4 


13.4 


77.1 


11.1 


-0.009 


0.158 


-0.121 


0.158 


-52.6 


25.0 


54.3 


32.0 
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Table 2 — Continued 



R 


Vo 


A t r 

AV 


O"0 


Acx 


h 3 


A 1 

Ah 3 


hi 


A 1 

Ah 4 


V ,c 


ATT" 

AV , C 




Aa , c 


(1) 


(2) 




(3) 




(4) 




(5) 




(6) 




(7) 




-1.044 


-44.9 


13.5 


74.9 


13.8 


0.079 


0.164 


-0.054 


0.164 


-34.6 


25.3 


65.0 


33.2 


-1.145 


-48.6 


14.9 


67.0 


11.8 


0.037 


0.202 


-0.132 


0.202 


-44.3 


27.8 


45.3 


35.4 


-1.247 


-54.2 


15.0 


73.8 


14.1 


0.064 


0.184 


-0.082 


0.184 


-46.0 


27.9 


59.0 


36.3 


-1.372 


-51.8 


23.3 


109.2 


23.9 


-0.044 


0.194 


-0.053 


0.194 


-60.1 


43.5 


95.0 


57.2 


-1.549 


-69.3 


17.1 


83.9 


12.4 


0.180 


0.186 


-0.154 


0.186 


-43.1 


32.2 


52.3 


40.5 


-1.778 


-52.1 


16.9 


75.6 


15.4 


-0.075 


0.204 


-0.092 


0.204 


-61.9 


31.7 


58.5 


40.9 


-2.052 


-56.6 


16.3 


81.1 


13.5 


0.077 


0.183 


-0.121 


0.183 


-45.8 


30.5 


57.1 


39.0 


-2.382 


-66.0 


18.0 


78.2 


19.2 


0.172 


0.209 


-0.039 


0.209 


-42.8 


34.0 


70.7 


44.4 


-2.765 


-54.6 


29.4 


102.9 


29.4 


0.048 


0.259 


-0.061 


0.259 


-46.1 


54.8 


87.5 


71.7 


-3.185 


-53.9 


21.7 


62.1 


17.2 


0.129 


0.317 


-0.131 


0.317 


-40.1 


40.6 


42.2 


51.5 


-3.673 


-70.0 


51.5 


167.0 


55.9 


-0.003 


0.280 


-0.033 


0.280 


-70.8 


96.0 


153.5 


127.6 


-4.206 


-46.9 


25.2 


30.0 


33.3 


0.039 


0.763 


0.047 


0.763 


-44.9 


47.0 


33.5 


65.3 


-4.722 


-67.4 


32.1 


6.5 


53.7 


-0.295 


4.505 


0.167 


4.505 


-70.7 


65.8 


9.1 


92.0 



Note. (1) Distance from center of M32 in arc seconds. (2) Rotation parameter V 
in km s -1 . (3) Dispersion parameter o"o in km s -1 . (4) LOSVD skewness parameter h 3 . 
(5) LOSVD kurtosis parameter J14. (6) Vo, c = Vo + v^3o"o^-3) an estimate of the true mean 
line-of-sight velocity. (7) a 0tC = a (l + y/Eh^), an estimate of the true line-of-sight velocity 
dispersion. 
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Table 3. M32 Kinematics as Derived via MPL 



R 


V 


AVo 


O"0 


Aa 


h 


A/i 3 


/l4 




V ,c 


AV ,c 


O"0,c 


Acr , c 


(1) 


(2) 




(3) 




(4) 




(5) 




(6) 




(7) 




4.559 


18.9 


12.8 


51.0 


10.4 


0.176 


0.140 


0.022 


0.033 


37.9 


19.7 


53.7 


12.1 


4.184 


70.2 


18.9 


97.2 


30.0 


-0.235 


0.155 


0.040 


0.079 


7.15 


51.4 


106.7 


31.4 


3.545 


47.3 


8.4 


47.9 


6.8 


0.054 


0.102 


0.008 


0.032 


56.9 


10.1 


48.8 


8.9 


3.068 


54.4 


8.1 


60.3 


12.7 


-0.054 


0.098 


0.094 


0.033 


49.5 


14.3 


74.2 


16.4 


2.632 


44.8 


6.7 


45.4 


7.0 


-0.008 


0.111 


0.028 


0.018 


35.4 


14.1 


48.6 


9.6 


2.247 


27.2 


9.8 


70.0 


10.1 


0.043 


0.094 


0.111 


0.060 


29.0 


12.9 


89.0 


14.4 


1.922 


37.6 


7.8 


62.7 


9.0 


-0.102 


0.080 


0.048 


0.031 


26.9 


11.4 


70.0 


10.1 


1.648 


71.8 


7.9 


64.6 


12.6 


-0.124 


0.076 


0.080 


0.027 


60.2 


16.7 


77.3 


16.9 


1.415 


46.4 


6.9 


54.8 


7.1 


-0.008 


0.079 


0.069 


0.023 


46.1 


10.0 


64.1 


9.5 


1.242 


44.7 


7.2 


64.4 


8.6 


-0.079 


0.091 


0.077 


0.026 


37.6 


10.7 


76.5 


11.6 


1.111 


49.1 


8.1 


87.5 


11.2 


0.070 


0.079 


0.037 


0.047 


62.7 


13.2 


95.5 


11.8 


1.009 


58.2 


7.2 


72.7 


11.2 


-0.144 


0.087 


0.104 


0.041 


37.1 


13.9 


91.1 


14.1 


0.908 


71.7 


4.6 


49.6 


6.5 


-0.034 


0.072 


0.036 


0.012 


69.3 


6.7 


54.0 


8.5 


0.806 


69.1 


4.9 


68.7 


8.3 


0.000 


0.070 


0.082 


0.023 


70.6 


8.2 


82.6 


8.7 


0.705 


49.5 


7.0 


67.4 


8.7 


0.138 


0.090 


0.071 


0.026 


74.2 


13.4 


79.1 


13.4 


0.654 


65.3 


5.3 


64.0 


7.1 


-0.004 


0.060 


0.041 


0.016 


63.5 


7.4 


70.4 


8.5 


0.604 


68.3 


5.2 


61.7 


7.7 


-0.157 


0.056 


0.051 


0.018 


52.0 


9.3 


69.4 


8.1 


0.553 


62.2 


5.4 


56.3 


6.4 


-0.034 


0.060 


0.060 


0.019 


61.5 


7.0 


64.6 


8.8 


0.502 


64.6 


4.4 


54.8 


6.7 


-0.022 


0.053 


0.070 


0.022 


63.9 


5.9 


64.2 


7.6 


0.451 


62.5 


5.3 


74.1 


7.0 


-0.071 


0.066 


0.054 


0.028 


53.6 


7.9 


83.9 


7.1 


0.401 


63.9 


5.0 


61.8 


7.6 


-0.040 


0.048 


0.089 


0.018 


61.8 


7.4 


75.3 


10.6 


0.350 


55.4 


4.3 


78.0 


7.5 


-0.040 


0.049 


0.062 


0.024 


51.1 


7.7 


89.8 


9.2 


0.299 


60.1 


4.5 


76.8 


7.0 


-0.178 


0.046 


0.053 


0.023 


37.4 


8.9 


86.8 


6.5 



Table 3 — Continued 



K 


T / 


AT/" 

AVo 


O"0 




h 3 


A U 

Ah 3 






Vo,c 


AT/ 

AVo, c 


O"0,c 




I 1 ) 


(2) 




(3) 




( 4 ) 




(5) 




(6) 




( 7 ) 




0.249 


62.5 


3.9 


66.5 


6.8 


-0.021 


0.039 


0.080 


0.029 


60.1 


5.2 


79.6 


7.2 


0.198 


55.2 


3.9 


82.4 


6.6 


-0.065 


0.042 


0.066 


0.022 


45.9 


6.7 


95.8 


6.2 


0.147 


43.3 


4.7 


101.2 


6.0 


-0.004 


0.048 


-0.022 


0.035 


42.5 


7.9 


95.8 


6.4 


0.096 


38.8 


5.1 


110.3 


6.3 


-0.009 


0.044 


0.032 


0.032 


41.0 


8.3 


118.9 


7.1 


0.046 


37.2 


5.0 


120.6 


7.0 


-0.028 


0.038 


0.118 


0.034 


31.3 


9.1 


115.6 


9.6 


-0.005 


-3.7 


4.7 


132.0 


6.9 


-0.051 


0.037 


0.139 


0.034 


-15.2 


10.7 


176.8 


11.5 


-0.056 


-24.5 


5.4 


130.0 


7.0 


-0.006 


0.037 


0.010 


0.034 


-25.7 


10.1 


133.2 


9.7 


-0.106 


-31.6 


4.7 


103.0 


6.0 


0.018 


0.039 


0.028 


0.030 


-28.1 


7.5 


110.2 


8.4 


-0.157 


-39.2 


4.9 


92.4 


6.8 


0.026 


0.042 


0.076 


0.028 


-35.9 


8.0 


109.7 


8.5 


-0.208 


-52.1 


4.8 


87.6 


6.5 


0.033 


0.041 


0.033 


0.031 


-47.5 


7.1 


94.8 


6.4 


-0.259 


-63.7 


4.9 


79.8 


7.0 


0.012 


0.042 


0.062 


0.030 


-62.1 


6.4 


91.9 


7.9 


-u.ouy 


PA 
-04. Z 


o.u 


pi p 

( .0 


P K 
0.0 


n nan 
u.uyu 


U.Uoy 


U.UOo 


u.uzz 


-00.0 


p 1 
0.0 


78 Q 


q n 
y.u 


-0.360 


-65.6 


5.8 


68.0 


10.4 


0.089 


0.049 


0.122 


0.030 


-56.0 


11.5 


88.3 


12.4 


-0.411 


-55.7 


5.5 


65.1 


7.8 


0.017 


0.043 


0.102 


0.022 


-54.3 


7.6 


81.5 


10.8 


-0.461 


-59.9 


4.9 


58.8 


6.3 


0.071 


0.055 


0.067 


0.025 


-53.7 


6.0 


68.4 


7.4 


-0.512 


-59.0 


5.4 


50.7 


6.4 


0.028 


0.046 


0.061 


0.018 


-58.4 


5.8 


58.3 


9.8 


-0.563 


-53.9 


6.5 


72.4 


8.5 


0.097 


0.061 


0.058 


0.025 


-40.3 


9.6 


82.7 


10.4 


-0.614 


-54.2 


7.4 


67.2 


12.0 


0.036 


0.064 


0.127 


0.040 


-51.1 


11.2 


88.0 


14.8 


-0.664 


-49.5 


6.0 


49.3 


7.6 


0.028 


0.057 


0.077 


0.030 


-48.4 


9.3 


58.5 


12.2 


-0.715 


-59.2 


7.1 


61.0 


9.4 


0.175 


0.061 


0.069 


0.028 


-42.6 


14.3 


71.3 


14.1 


-0.816 


-50.8 


6.8 


55.6 


9.5 


0.042 


0.055 


0.113 


0.029 


-49.6 


10.5 


71.0 


13.8 


-0.918 


-43.3 


6.9 


58.9 


10.1 


0.050 


0.080 


0.101 


0.031 


-36.2 


10.3 


73.5 


15.2 


-1.019 


-50.7 


5.6 


60.6 


6.8 


-0.004 


0.065 


0.076 


0.041 


-50.7 


8.2 


72.0 


9.1 


-1.121 


-56.2 


6.7 


53.3 


8.1 


0.080 


0.073 


0.066 


0.027 


-49.7 


14.2 


61.9 


12.6 
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Table 3 — Continued 



R 


V 


AV 


00 


A(T 


h 3 


A/i 3 


hi 


Ah 4 


V ,c 


AVb )C 


O"0,c 


Ao- , c 


(1) 


(2) 




(3) 




(4) 




(5) 




(6) 




(7) 




-1.252 


-58.6 


7.3 


60.2 


8.5 


0.092 


0.072 


0.066 


0.022 


-47.2 


10.8 


69.9 


12.4 


-1.425 


-55.7 


9.5 


88.1 


18.1 


0.035 


0.077 


0.182 


0.061 


-48.3 


17.0 


127.4 


21.6 


-1.658 


-64.8 


6.4 


51.2 


7.4 


0.037 


0.058 


0.056 


0.019 


-60.5 


9.0 


58.2 


9.4 


-1.932 


-56.9 


8.1 


58.0 


8.1 


0.009 


0.094 


0.076 


0.041 


-54.9 


9.8 


68.8 


9.9 


-2.257 


-75.5 


9.5 


67.2 


13.5 


0.164 


0.096 


0.084 


0.041 


-43.6 


17.1 


81.1 


20.4 


-2.642 


-57.6 


10.3 


62.1 


11.9 


0.037 


0.105 


0.068 


0.036 


-47.4 


15.9 


72.5 


16.9 


-3.078 


-56.0 


9.9 


53.8 


11.1 


0.268 


0.106 


0.088 


0.043 


-28.7 


17.7 


65.4 


15.9 


-3.555 


-41.6 


27.4 


75.6 


50.9 


0.070 


0.283 


0.404 


0.302 


-28.8 


34.7 


150.3 


30.9 


-4.194 


-9.2 


15.6 


56.4 


14.9 


-0.144 


0.105 


0.021 


0.037 


-22.5 


34.1 


59.3 


20.1 


-4.569 


-39.7 


10.8 


45.3 


9.2 


0.240 


0.136 


0.042 


0.036 


-19.7 


17.4 


49.9 


13.2 
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Fig. 1.— 

Top: STIS spectra of M32 at five positions along the slit. Solid curves are convolutions 
of the MPL-derived broadening functions N(V) (Fig. with the stellar template. Bottom: 
Spectrum of HR7615, a KOIII giant, the template spectrum. The vertical scale of the 
template spectrum is compressed with respect to that of the M32 spectra. 

Fig. 2.- 

Broadening functions recovered from the central spectrum of M32 using the two 
spectral deconvolution algorithms discussed in the text. The degree of smoothing increases 
downward. Left column: FCQ. (a) W = 2.4; (b) W = 1.6; (c) W = 1.2; (d) W = 0.8. Right 
column: MPL. (e) a = 10 3 ; (f) a = 10 5 ; (g) a = 10 7 ; (h) a = 10 9 . The MPL estimates 
tend toward a Gaussian for large a while the FCQ estimates become increasingly distorted 
as the smoothing is increased. This is the source of the greater bias in the FCQ estimates 
(although in practice smoothing parameters as small as W = 0.8 would never be used). 

Fig. 3.- 

The mean value of computed by the FCQ algorithm between 0.5" and 2.0", as a 
function of smoothing parameter W. Large values of W correspond to small degrees of 
smoothing and hence to less biased estimates. These h 4 values are based on data that were 
heavily binned in radius in order to increase the S/N as much as possible. 

Fig. 4.— 

Line-of-sight velocity distributions derived from the STIS M32 spectra using the MPL 
deconvolution algorithm. Note the sudden increase in the width of the broadening functions 
inside of ~ 0.2". The LOSVDs are roughly antisymmetric about the center of M32, as 
expected for a relaxed system; the right column shows N(V) averaged over the left and right 
sides, ±[N(V,R) + N(-V, -R)]. The central LOSVD exhibits strong non-Gaussian wings, 
a likely consequence of high-velocity stars near the central black hole. The broadening 
functions at larger radii exhibit asymmetries suggestive of a second kinematic subcomponent 
which rotates with a velocity closer to the systemic velocity of M32. 

Fig. 5.- 

STIS rotation curve for M32, derived from LOSVDs obtained using the MPL spectral 
deconvolution algorithm. Upper panel: filled circles: Vo, the parameter that measures 
the velocity shift of the Gaussian function that multiplies the Gauss-Hermite series. 
Open circles: V + \/3a h 3 , an estimate of the true mean line-of-sight velocity. Lower 
panel: the Gauss-Hermite parameter h% that measures asymmetries in the LOSVDs. The 
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mean velocity is smaller than |Vo| due to the nonzero value of I13, which in turn reflects 
asymmetries in the LOSVD's (Fig. 

Fig. 6.— 

STIS velocity dispersion profile for M32, derived from LOSVDs obtained using the 
MPL spectral deconvolution algorithm. Upper panel: filled circles: (Jq, the parameter that 
measures the dispersion of the Gaussian function that multiplies the Gauss-Hermite series. 
Open circles: o"o(l + Vhl), an estimate of the true rms line-of-sight velocity. Lower panel: 
the Gauss-Hermite parameter /i 4 that measures the amplitude of symmetric non-Gaussian 
distortions in the LOSVD. The velocity dispersion is generally greater than cr due to the 
nonzero values of h^. This difference is substantial in the inner ~ 0.2" due to the strongly 
non-Gaussian wings of the central LOSVDs (Fig. |]). 

Fig. 7.- 

Comparison of Vq and <7o derived from the M32 STIS data (filled circles) with earlier 
ground-based determinations. Squares: WHT measurements from van der Marel et al. 
(1994a). Triangles: CFHT measurements from Bender et al. (1996). 

Fig. 8.— 

Comparison of Vq and do derived from the M32 STIS data (filled circles) with FOS 
data of van der Marel et al. (1997) (open squares). 

Fig. 9.— 

Reduced \ 2 contours describing the fit of the axisymmetric models described in the 
text to the observed, mean square line-of-sight velocity at points within the inner arc second 
of M32. The contour spacing is 0.5 and the innermost contour is at 0.75. The plus symbol 
marks the best-fit model. 

Fig. 10.— 

Predicted kinematical profiles for three axisymmetric models with different black hole 
masses, (a) M h = 2.0 x 1O 6 .M ; (b) M h = 3.0 x 1O 6 .M ; (c) M h = 4.0 x 10 6 M G . Thin 
curves show the models as observed with infinite resolution; heavy curves are the models 
after convolution with the STIS PSF; open circles are the data points. For each M^, the 
mass to light ratio M/L and rotational parameter k have been adjusted to optimize the fit. 

Fig. 11 — 

The central M32 LOSVD (heavy line), symmetrized about V = 0, compared to the 
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LOSVD's predicted by spherical nonrotating models with three black holes masses. Thin 
line: M h = 2.5 x 1O 6 .M ; dashed line: M h = 5.0 x 10 6 M Q ; thin line: M h = 10. x 10 6 M Q . 

Fig. Al - 

Estimated Gauss-Hermite parameters a and /i 4 , derived by fitting the N(V) of equation 
(A6) to the assumed form (A5), with j max = 4. Filled circles indicate the input values of a 
and /i 4 ; these values are recovered only when the input N(V) has h 6 = 0. 

Fig. Bl — 

Recovery of a via FCQ. 

Fig. B2 — 

Recovery of fo 4 via FCQ. Circles: do = 40 km s~ x . Triangles: <7o = 60 km s _1 . 
Diamonds: a = 100 km s -1 . There is a significant negative bias in the recovered values of 
/i 4 when the velocity dispersion is less than about 100 km s -1 . 

Fig. B3 - 

Mean estimates of N(V) averaged over 100 random realizations of the observed 
spectrum, for S/N = {5,20, 100}. The input N(V) (Equation B2) is shown by the heavy 
curves. 

Fig. B4 — 

MISE and ISB of estimates of N(V) obtained from the two deconvolution algorithms. 
The input LOSVD was a Lorentzian (Equation B2) with a = 108km s" 1 and /i 4 = 0.15. 
Both MISE and ISB have been normalized as described in the text. Solid lines: MPL 
algorithm. Open circles: FCQ algorithm, using a fixed smoothing parameter W = 1.3. 
Filled circles: FCQ algorithm, using the value W opt that minimizes the MISE of the 
estimated N(V). 

Fig. B5 - 

MSE and bias in estimates of h 4 obtained from the two deconvolution algorithms. The 
input N(V) was a Lorentzian (Equation B2) with er = 108km s" 1 and h 4 = 0.15. Solid 
lines: MPL algorithm. Open circles: FCQ algorithm, with W = 1.2. Filled circles: FCQ 
algorithm, using the value W opt that minimizes the MISE of the estimate h A . W opt is plotted 
vs. S/N in the bottom panel. 
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Fig. I.- 




Fig. 2.- 
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Fig. 5.- 
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Fig. 9.- 
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Fig. 10.- 
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Fig. 11.- 
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Fig. B3 
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